DCP-10-02 



Implications of Yukawa Textures in the decay H + — > W + 'y within the 2HDM-III 

J. E. Barradas-Guevara, 1 '!*] F. Cazarez-Bush, 2, [ll A. Cordero-Cid, 3 '^ 0. 

Felix-Beltran, 3 ^ J. Hernandez-Sanchez, 3 ' 4 'Q and R. Noriega-Papaqui 5 ' 4 '[** 

J Fac. de Cs. Fisico- Matemdticas, BUAP. Apdo. Postal 1364, C.P. 72000 Puebla, Pue., Mexico. 
2 Instituto de Matemdticas, Universidad Nacional Autonoma de Mexico, 
Cd. Universitaria, Mexico, D.F. 01000, Mexico. 
3 Facultad de Ciencias de la Electronica, Benemerita Universidad Autonoma de Puebla, 
Apdo. Postal 542, C. P. 72570 Puebla, Pue., Mexico. 
^Dual C-P Institute of High Energy Physics, Puebla, Pue., Mexico. 
5 Centro de Investigacion en Matemdticas, 
Universidad Autonoma del Estado de Hidalgo, 
Carr. Pachuca-Tulancingo Km. 4-5, CP. 4^184, Pachuca, Hgo., Mexico. 
(Dated: February 26, 2010) 

We discuss the implications of assuming a four-zero Yukawa texture for the charged Higgs 
decay H + — > W +r y at one-loop level, within the context of the general 2-Higgs Doublet Model 
of Type III. We begin by presenting a detailed analysis of the charged Higgs boson couplings 
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I. INTRODUCTION 

Detecting a charged Higgs boson during the imminent Large Hadron Collider (LHC) experi- 
mental running would constitute a clear evidence of physics beyond the Standard Model (SM) 
Charged Higgs bosons appear in many well motivated extensions of the SM, whose phenomenology 



has been widely studied over the years 
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ar, 2-Higgs Doublet Models (2HDMs), in 
can be considered as a prototype of a 



iC, it will be a future International Linear Collider (ILC) [9] or Compact Linear Collider 
ip| which will have the definite word about exactly which mechanism of mass generation 



-4|. In particu 
both Supersymmetry (SUSY) and non-SUSY versions 
Higgs sector that includes a charged Higgs boson (H^). It is expected that the LHC will allow us 
to test the mechanism of Electro- Weak Symmetry Breaking (EWSB) and, in particular, to probe 
the properties of charged Higgs bosons, which represent a unique probe of a weakly-interacting 
theory, as is the case of the Minimal Supersymmetric Standard Model (MSSM) [5] and general 
2HDMs of Type I, II, III and IV (2HDM-I, 2HDM-II, 2HDM-III and 2HDM-IV) [7], or whether 
strongly-interacting scenarios are instead realized, like in the old Technicolor models or the ones 
discussed more recently jj]. Ultimately, while many analyses in this direction can be carried out 
at the L 
(CLIC) 

and which realization of it occurs in Nature. 

The 2HDM-II has been quite attractive to date, in part because it coincides with the Higgs 
sector of the MSSM, wherein each Higgs doublet couples to the u- or ci-type fermions separately 1 . 
However, this is only valid at tree-level [12 ]. When radiative effects are included, it turns out 
that the MSSM Higgs sector corresponds to the most general version of the 2HDM, namely the 
2HDM-III, whereby both Higgs fields couple to both quarks and leptons. Thus, we can consider 
the 2HDM-III as a generic description of physics at a higher scale (of order TeV or maybe even 
higher), whose low energy imprints are reflected in the Yukawa coupling structure. With this idea 
in mind, a detailed study of the 2HDM-III Yukawa Lagrangian was presented in Refs.fl^l. [ul ]. 
under the assumption of a specific texture pattern which generalizes the original model of 
Ref. [16]. Phenomenological implications of this model for the neutral Higgs sector, including 
Lepton Flavour Violation (LFV) and/or Flavour Changing Neutral Currents (FCNCs) have been 
presented in a previous work [17J] . The extension of such an approach to investigate charged Higgs 
boson phenomenology was conducted in Ref. [14], which discussed the implications of this Yukawa 
texture for the charged Higgs boson properties (masses and couplings) and the resulting pattern 



1 Notice that there exist significant differences between the 2HDM-II and MSSM though, when it comes to their 
mass/coupiing configurations and possibie Higgs signais [Tlj ] . 
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of charged Higgs boson decays and main production reactions at the LHC. 



Decays of charged I 
modes W ± -/,W ± Z° [18, 



i ggs bosons have been studied in the literature, including the radiative 
1 1 Of ] . mostly within the context of the 2HDM-II or its SUSY incarnation 



(i.e., the MSSM), but also by using an effective Lagrangian extension of the 2HDM 20|. More 
recently, within an extension of the MSSM with one Complex Higgs Triplet (MSSM+1CHT) [21 . 



221 ] . While the decay mode H + — > is forbidden at the tree level due to electromagnetic 

gauge invariance, the H + — > W + Z decay can be induced at this order in models including Higgs 



triplets or more complicated representations 



211 ] . In spite of their suppressed branching ratios, 



these decay modes are very interesting, due to the fact that their experimental study may provide 
important information concerning the underlying structure of the gauge and scalar sectors; these 
channels have a clear signature and might be at the reach of current and future particle colliders. 



Charged Higgs boson production at hadron colliders was studied long ago [2i 



ana, more 



recently, systematic calculations of production processes at the LHC have been presented 



25]. 



Current bounds on the mass of a charged Higgs boson have been obtained at Tevatron, by 



studying the top decay t — > bH + , which already eliminates large regions of t 



ie_parameter space 



27 



28j, rather model 



261 ]. whereas LEP2 bounds imply that, approximately, m H + > 80 GeV 
independently. Concerning theoretical limits, tree-level unitarity bounds on the 2HDM Higgs 
masses have been studied in generic 2HDMs and in particular an upper limit for the charged Higgs 



22]- 



mass of 800 GeV or so can be obtained, according to the results of Ref. 

In this paper we extend previous studies, to include the decay H + — > W + j. This paper 
is organized as follows. In section II, we discuss the Higgs- Yukawa sector of the 2HDM-III, in 
particular, we derive the expressions for the charged Higgs boson couplings to heavy fermions. 
Then, in section III, we review the expressions for the decays H + — > VF + 7 at one-loop level 
and numerical results are presented for some 2HDM-III scenarios, defined for phenomenological 
purposes. Actual LHC event rates for the main production mechanisms at the LHC are given 
in section IV. These include the s-channel production of charged Higgs bosons through cb(cb)- 
fusion [3^] and the multi-body more qq, gg — > tbH~+ c.c. (charge conjugated). These mechanisms 
depend crucially on the parameters of the underlying model and large deviations could be expected 
in the 2HDM-III with respect to the 2HDM-II. Finally, we summarize our results and present the 
conclusions in section V. 
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II. THE CHARGED HIGGS BOSON LAGRANGIAN AND THE FERMIONIC 

COUPLINGS 



We shall follow Refs. 1_3|, [lTj] , where a specific four-zero texture has been implemented for the 
Yukawa matrices within the 2HDM-III. This allows one to express the couplings of the neutral and 
charged Higgs bosons in terms of the fermion masses, Cabibbo-Kobayashi-Maskawa (CKM) mixing 
angles and certain dimensionless parameters, which are to be bounded by current experimental 
constraints. Thus, in order to derive the interactions of the charged Higgs boson, the Yukawa 
Lagrangian is written as follows: 



C Y = Y?Q L $ 1UR + Y 2 u Q L $ 2 u R + YfQL^dR + Y 2 d Q L $ 2 d R , 



(1) 



where = (<l ) i2i ( t ) i2) T refer to the two Higgs doublets, $>i j2 = ia 2 ^i 2 , Ql denotes the left- 
handed fermion doublet, u R and d R are the right-handed fermions singlets and, finally, Y^ 2 denote 
the (3 x 3) Yukawa matrices. Similarly, one can write the corresponding Lagrangian for leptons. 

After spontaneous EWSB and including the diagonalizing matrices for quarks and Higgs 
bosons 2 , the interactions of the charged Higgs boson H + with quark pairs have the following 
form: 



H+ 



2V2M- 



w l=1 1 



CKM )il 



tan (3 m<2; — sec j3 



+ 



cot j3 m Ui 5u - esc j3 



( V2M] 



V 9 



+(Vckm)u 



tan (3 Sij — sec f3 




(2) 



cot P m Ui Su - esc P ^^Mk 

where Vckm denotes the mixing matrices of the quark sector (and similarly for the leptons). The 
term proportional to Sy corresponds to the contribution that would arise within the 2HDM-II, 
while the terms proportional to Y 2 d and YJ" denote the new contributions from the 2HDM-III. 
These contributions, depend on the rotated matrices: Yn = OT P q Yn Pq O q (n = 1 when q = u, 
and n = 2 when q = d ), where O q is the diagonalizing matrix, while P q includes the phases of 
the Yukawa matrix. In order to evaluate YS we shall consider that all Yukawa matrices have the 



2 The details of both diagonalizations are presented in Ref. 



four-Hermitic-texture form [151 ] . and the quark masses have the same form, which are given by: 



( o c q o N 



M q 



\ 



C* Bg Bg 

B* A q ) 



(q = u,d). 



(3) 



To diagonalize these matrices, we use the matrices O q and P q , in the following way [la ]: 



m = T q PgMlP\Og. 



(4) 



Then, one can derive a better approximation for the product Oj P q Y q Pq O q , expressing the rotated 
matrix iv, in the form 



Y q 



m q m q 



[Xrilii 



mfrrij 



[Xrilij & 



(5) 



In order to perform our phenomenological study, we find it convenient to rewrite the Lagrangian 
given in Eq. ([2]) in terms of the coefficients [xnjyj a s follows: 

3 



dklH (FcKM) 



n . sec/3 a 

tan p m dl dy j=- ^/m dl m dj Xij 



cot P m Ui 5 a - < ^ 7 =- ^Jm Ui m Ul Xu 



+ 

+{Vckm)u 



V2 



SGC j3 

tan/3m di <% j=- y/m dl m dj Xij 



V2 

(VcKM.)lj 
d 



V2 



1 



cot P m Ui S u - ^i- y/m Ui m Ul Xa 



(VcKM)ij 7 > dj H 



(6) 



where we have redefined [X\]ij = Xij an d = X%- Then, from Eq. ([6]), the couplings uidjH' 

and UidjH~ are given by: 



QH+UiC 



w 



-(Sij + Pij%), g H - ud 



'9 



2V2M W 



(Sij - P«75), 



(7) 



where Sij and Py are defined as: 



Sa = E \^ 



i=i 
+ 

3 



'CKMjiZ 



sec P d 

tan/^m rf; d tj j=- ^m dl m dj Xij 



V2 



esc P 



Pa = £ m 



1=1 



CKM)il 



V2 

tan Pm dl Sy 



cot /3 m Ui 5u j=- s /m Ui m Ul Xu 

sec P / 



V2 



i c CSC (3 

cot P m Ui Su - —j=- yjm ui m n Xu 



(Vckm)o| 
y/m dl m d . xtj 
(Vckm)«j|- 



(8) 



As it was discussed in Ref. [13j, most low-energy processes imply weak bounds on the coefficients 
xfj, which turn out to be of OCl). However, some important constraints on tan/3 have started to 
appear, based on l?-physics [311 ]. In order to discuss these results we find convenient to generalize 
the notation of Ref. [32] and define the couplings UidjH + and UidjH~ in terms of the matrices 
Xij, Yij and (for leptons). In our case these matrices are given by: 

sec/3 [mi 



X, 



tan /3 Sij 
cot /3 8u - 



V2 

esc /3 



m d . 



3 xt 



in 



u i ~u 



Xu 



(9) 



where Xu and Yu are related with Sij and Py defined in the Eq. (jSJ) as follows: 



Sij = ^(Vckm)^ m di Xij + m Ui Yu(V CKM )ij, 
i=i 

3 

Pij = y^ J { v CKu)u m dl Xij - m Ui Yu(y C KM)ij- 



(10) 



i=i 



The 33 elements of these matrices reduce to the expressions for the parameters X,Y,Z 
X33, Y33, Z33) used in Ref.[32j]. Based on the analysis of B — > 



32 



33j , it is claimed that X < 20 



and Y < 1.7 for mjy+ > 250 GeV, while for a lighter charged Higgs boson mass, mjj+ ~ 180 GeV, 
one gets (X, Y) < (18,0.5). In recently work we get the values of (X,Y) as a function of tan/3 
within our model. Thus, we find the bounds: Ix^l ~ 1 f° r 0-1 < tan/3 < 70 Although in 

our model there are additional contributions (for instance from c-quarks, which are proportional 
to X23), they are not relevant because the Wilson coefficients in the analysis of B — >■ Xgj are 
functions of m^/My^ or w? c /m 2 H+ [3J, that is, negligible when compared to the leading X33 effects, 
whose Wilson coefficients depend on m^/My^ or m1/m 2 H+ . Other constraints on the charged Higgs 
mass and tan/3, can be obtained the anomalous magnetic momentum of the muon based on Aa^, 



35|, 



36|. For instance, 



the p parameter, as well as B-decays into the tau lepton, can be obtained 
as can be read from Ref. [37], one has that the decay B — > tu, implies a constraint such that for 
m,H+ = 200 (300) GeV, values of tan/3 less than about 30 (50) are still allowed, within MSSM 
or THDM-II: However, these constraints can only be taken as estimates, as it is likely that they 
would be modified for THDM-III. In summary, we find that low energy constraints still allow to 
have xt; = O(l) 3 . On the other hand, the condition rg+ < i in the frame of the 2HDM-II implies 
which leads to 0.3 % tan/3 < 130. However, in the 2HDM-III we have that 



H+ 



3G F mf 



l H+ AV2-K tan /3 2 



A more detailed analysis that includes the most recent data is underway [3a . I3S 
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H+ 



AV2-K tan (3 2 



when IX33 
2HDM-II |7, 



V2c 



we have checked numerically that this leads to 0.08 < tan j3 < 200 



1 and 0.3 < tan j3 < 130 as long as IX33I ~~ * recovering the result for the case of the 



401 ]. In this sense, if we consider the constraints imposed by the perturbativity bound, 



a portion of the low tan j3 appearing in some graphs would be excluded. However, we have decided 
to keep that range both to show the behaviour of the quantities of interest, and also to keep in 
mind that such criteria (perturbativity) should be taken as an order of magnitude constraint. 



III. DECAY H+ -> W+-/ 

In the calculation of the width for decay H + — > W + 7, we employe a nonlinear i?£-gauge, which 
leads to considerably simplifications due to the fact that some unphysical vertices are removed from 
the interaction Lagrangian [3]. We have shown that such a gauge not only reduces considerably 
the number of Feynman diagrams but also renders manifestly gauge-invariant and ultraviolet-finite 
amplitudes. Apart from emphasizing the advantages of using the nonlinear i?g-gauge, we will 
analyze the H + — > W + j decay in some scenarios which are still consistent with the most recent 
bounds discussed in the previous section. 



A. The gauge fixing procedure 

We now would like to comment on the gauge-fixing procedure which was used to simplify our 
calculation. To this end we introduce the following gauge-fixing functions [411 ] : 

/ + = ( D l + W + » - itmwG+r, (11) 

f = - imzGz, (12) 

f A = d^, (13) 

with the electromagnetic covariant derivative and £ the gauge parameter. Note that / + is 
nonlinear and transforms covariantly under the electromagnetic gauge group. This gauge- fixing 
procedure is suited to remove the unphysical vertices WGwl and WG\yZ, which arise in the Higgs 
kinetic-energy sector, and also modifies the Yang-Mills sector. One important result is that the 
expression for the WW'j vertex satisfies a QED-like Ward identity, which turns out to be very 
useful in loop calculations. In particular, in the Feynman-t'Hooft gauge the Feynman rule for the 
W+(q)W~ (p)A^(k) coupling can be written as 

rj^T 7 = ~ ie (g^(p ~ k + q) P + 9»M - v)» + sv^ -v- <?W , (14) 
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where all the momenta are incoming. It is easy to see that this expression satisfies a QED-like 
Ward identity. The remaining Feynman rules necessary for our calculation do not depend on the 
gauge fixing procedure. As far as the Yukawa couplings are concerned, we will concentrate on the 
type-Ill 2HDM. 



B. Decays of the charged Higgs boson at tree level 



The expressions for the charged Higgs boson decay widths H 

3<? 2 



Uidj are of the form: 



r(ir 



mdj) 



32ttM?, 



W m H+ 



\l/2/ 2 2 2 \ 

-A ' (m H+ ,m Ui ,m d ) 



in 



H+ 



711 „ 



(Sfj + Pi 



, "'« "'-/, (-^r, 



(15) 



where A is the usual kinematic factor A(a, b, c) = (a — b — c) 2 — 46c. When we replace Xud — > 0, 
the formulae of the decays width become those of the 2HDM-II: see, e.g., Ref. Furthermore, 



the expressions for the charged 
as in the 2HDM-II. Recently 



3iggs boson decay widths of the bosonic modes remain the same 
141 ]. we have studied that the effect of the modified Higgs couplings 
typical of the 2HDM-III shows up clearly in the pattern of charged Higgs boson decays, which can 
be very different from the 2HDM-II case and thus enrich the possibilities to search for states 
at current (Tevatron) and future (LHC, ILC/CLIC) machines. 



C. Decay H + — » W + "f induced at one-loop level 

In the nonlinear i?£-gauge, the decay H + — > W + j receives contributions from the Feynman 
diagrams shown in Fig. [TJ As far as the fermionic sector is concerned, the main contribution 
comes from the third-generation quarks (t,b), which induce three diagrams. However, we also 
considered the quarks contribution (c, b) and (t, s) because their coupling with the charged Higgs 
can be important [bj ]. which induce six more diagram. Whereas in the bosonic sector there are 
contributions from the (H^, (fp) and (VF 1 * 1 , <jP) pairs, with cfp = h° or H°. We would like to 
emphasize that the nonlinear i?g-gauge considerably simplifies the calculation of the decay H + — > 
V + j. First of all, the removal of the unphysical vertex W^G^yj, and the tadpole graphs vanish 
Ijj]. Apart from these simplifications, we will be able to group the Feynman diagrams into subsets 
which separately yield a manifestly gauge-invariant amplitude free of ultraviolet singularities. Once 
the amplitude for every Feynman diagram shown in Fig. [1] was written down, the Passarino- 
Veltman method |42] was applied to express it in terms of scalar integrals, which are suitable for 



FIG. 1: Feynman diagrams contributing to the H + — > W + ^f decay in the nonlinear i?^-gauge. <fr° stands for 



h° and H°, and X+ for H+ and G^. 



n 

numerical evaluation [43]. The full amplitude for the decay H + (p) — > W + (q)j(k) can be written 



MS 



3/2 

M{H+ -> W+ 1 ) = 2 J^ mwg2 (G(q-k 9ta , - Kq^) -iHe^k a q^ e^k)e u (q), (16) 

which is manifestly gauge invariant. The H function only receives contributions from the quarks: 

i 3 

H = 7T ( 3 ^ B {H± t W,dj,u.i) + A (H ± , W) y4j C {dj,Ui) ~ 2 Kj C ( Ul ,dj) ) ) > ( 17 ) 

where 



(18) 



Kj = ^(Pij + Sij), (19) 

Aj = ^m d< (^-^), (20) 

here, we have introduced the shorthand notation A( a fe ) = m^ — m 2 , AB( abcd ^ = Bo (m 2 , m 2 , m 2 ,) — 
Bq (m 2 ,m 2 ,m 2 d ^ , 4 and C( m = Co {jn 2 H ± , m^y, 0, mj, m 2 , m 2 ) , where Bq and Co stand for 
Passarino-Veltman scalar integrals. From the AJ3uf, }Cj( n definition, it is clear that the H func- 
tion is ultraviolet finite. Again, when we replace Xud 0, the formulae of the decays width 



become those of the 2HDM-II 



19]. As for C, it can be written as 



G = WZTi — X2 ( Gu ^ + ( G E±<p° + G G ± + G w ±^o) ) , (21) 

tm H± l\ {H± w) \ ^ 0= ^o jh o w 1 



When a = 0, m a must be substituted by 0. 
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where Gab stands for the contribution of the (A, B) pair: 

Gutdj = y^m^ ± ( (3 A (M - )d .) + m^] Xjj + 3ajjA {H ± w) jAB {H ± Wd . Ui) 

ij 

+ ( 2 «4.Ay - Xfj A (H ± jW) ) C( u . jd .^\ , (22) 

where 

ff u = \ ( m «, ( p v + s a) + - Sun ' ( 23 ) 

When one takes XijS - ^ 0, the formulas for the function G Ui ^. reduce to the 2HDM-II case, see, e. 



Gn+h = K H±y^(H±,h )^(H ± ,W) SB (0,W,H ± ,h ) ~ ( 2m -//±^(//±,h°) + m hO m w) d B (H±,W,H±,h ) 

- m 2 H± A {H ± tW) (l + 2m 2 H ±C (H ± >h oy) J , (24) 



G G± h° = fc° ( ( m VK A (/i°, W) + m H± (3m^ - 2m 2 h o)) &B(ff±,w,ft°,W) + ^{H±, h°) w)$ B < 



W,h°,W) ' ^(H±, ft. ) ^(/i , VI/JO-DCO.H/^ ,^) 



+ m^iA^i.iv) (l + 2m^C (H , ift o)) J, (25) 

and 

G w ± H o = \ w ±h°( (m^A (h o jW r) +m^± (7m^ - 2m^„ - 4m#±)) 5B {H ±^ H ± h0) 

+ A(H±,W)A(h°,W)$ B (0,W,H±,h o ) 

with 



+ rn H± A {H ± >w) (1-2 (2A {H±>w) -m w )C {m ^ m 2^ ) j J, (26) 



<\ff±/i° = c a -/3 2 s a _,3 m^± + /zf 2 + ( Sasptp - ) m 2 h0 , (27) 



s 2(a-P)A(H±,h°) 



(28) 



and 



m 2 w s 2{a -0) 

*W±h° = ^ ' ^ ' 

Finally, the contribution of the heaviest CP-even scalar boson H° is obtained of the lightest one, 
once the substitutions m h o — > m H o and a — > a — ir/2 are done. 
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It is also evident that the partial amplitudes induced by the pairs (t, b), (t, s), (c, b), (0°, W^), 
H^), and Gyy) are gauge invariant and ultraviolet finite on their own. Again, this is to 
be contrasted with the situation arising in the linear i?g-gauge, where showing gauge invariance 
is somewhat cumbersome, and the cancellation of ultraviolet divergences in the bosonic sector is 
achieved only after adding up all the Feynman diagrams 19]. 



Once Eq. (|16p is squared and the spins of the final particles are summed over, the decay width 
can be written as 



r(ii+ -+ w^) = f 2 f {H± 2 w \ (Id 2 + \H?) • (30) 



We will evaluate this decay width for some values of the parameters of the model. 



D. Numerical results and discussion 



Let us now discuss the decay modes of the charged Higgs boson within our model. Hereafter, we 
shall refer to two benchmark scenarios, namely, (i) Scenario A: xfj = ~~ 1> xfj = ~ lj (h) Scenario 
B: xfj = 1) Xij = 1- We have performed the numerical analysis of charged Higgs boson decays by 
taking tan/3 = 0.3, 0.5, 1, 10 and varying the charged Higgs boson mass within the interval 100 
GeV < m H ± < 800 GeV, further fixing m h o = 120 GeV, m A o = 300 GeV and the mixing angle 
at a = 7r/2. Then the results for the Branching Ratios (BRs) are shown in Fig. [2j and have the 
following characteristics. 

Scenario A. In Fig. [2]^ a) we present the BRs for the channels H + — > W + j, H + — > tb, cb, ts, t + v t , 
W + h°, W + A° as a function of m H +, for tan/3 = 0.3 and fixing m h o = 120 GeV, m A o = 300 GeV 
and the mixing angle a = n/2. When mjj+ < 175 GeV , we can see that the dominant decay of 
the charged Higgs boson is via the mode ti>, with BR(^ + — > rv) ~ 1, and the same importance 
we have the mode cb. For the case 170 GeV < < 180 GeV the mode ts is relevant. In this 

range, the mode VF + 7 is competitive and becomes of order 10 _1 , which it is also very different from 
the 2HDM-II and becomes an interesting phenomenological consequence of the 2HDM-III. We can 
also observe that, for nrijj+ > 180 GeV, the decay mode tb is dominant (as in the 2HDM-II). In 
this range mass of the charged Higgs, the channel H + — > W + 'y has a BR ~ 10 -4 to ~ 10~ 6 . Now, 
from Fig. E)^b), where tan/3 = 0.5, we find that the dominant decay mode is into t + v t for the 
range m H + < 175 GeV, again for 175 GeV < m H + < 180 GeV the mode ts is the leading one, 
but for 180 GeV< m H +, the decay channel tb becomes relevant, whereas for the range 170 GeV 
< m H + < 180 GeV the mode W + j has a BR of order 10~ 2 . Again, when m H + > 180 GeV the 
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FIG. 2: The figure shows the BRs of the H + decaying into the principal modes in Scenario A, taking 
Xf. = -1, xf. = -1, m h0 = 120 GeV, m A o = 120 GeV and a = tt/2 for: (a) tan/3 = 0.3, (b) tan/3 = 0.5, 
(c) tan/3 = 1, (d) tan/3 = 10. The lines in each graph correspond to: (Wj line) BR(iJ + — > VF + 7), (i& 
line) BR(iJ+ -> f£), (c6 line) BR(#+ -> c6), (is line) BR(#+ ts), (rv line) BR(£f+ r+^ r ), (Wh line) 
BR(fl+ -> W+fe ), (WA line) BR(iJ+ -> VK+A°). 

BR{H + —7- VF + 7) is of order 10 -4 to 10~ 7 . Then, see Fig. Wis)-, f° r the case with tan/3 = 1 one 
gets that BR(H + — > t + u t ) 1 when m#+ < 175 GeV, while in the range 175 GeV < m H + < 180 
GeV the mode ts is relevant and the BR(H + — > W + j) even could be important with an order 
10 -3 . However for 180 GeV< m H +, the dominant decay of the charged Higgs boson is the mode 
tb. For tan f3 = 10, we show in plot Fig. [2jd) that the dominant decay of the charged Higgs boson 
is the mode t + v t , when m H + < 180 GeV, but that, for 180 GeV< m H + < 250 GeV, the decay 
channel tb becomes the leading one, whereas for the range 250 GeV < m H +, the mode W + h° is 
relevant. In this case the BR of the decay channel W +r y is unimportant. 
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FIG. 3: The same as in Fig. H but taking \% = 1, X% = 1 (Scenario B). 

Scenario B. In Fig. [3](a) we present the BRs for the channels , H + — > W + ^,tb, cb, ts, t + v t , 
W + h°, W + A° as a function of m H +, for tan/3 = 0.3 and fixing m h o = 120 GeV, m A o = 300 GeV 
and the mixing angle a = tt/2. When < 175 GeV , we can see that the dominant decay of 

the charged Higgs boson is via the mode t + v t , with BR(.£f 4 + — > t + v t ) ~ 1. For the case 175 GeV 
< m H + < 180 GeV the mode ts is dominant, whereas for the range 170 GeV < m H + < 180 GeV 
the mode VF + 7 has a BR of order 10~ 2 . We can also observe that, for > 180 GeV, the decay 
mode tb is dominant (as in the 2HDM-II). Now, from Fig. 0(b), where tan/3 = 0.5, we find that the 
dominant decay mode is into t + v t for the range m H + < 175 GeV, again for 175 GeV < m H + < 180 
GeV the mode ts is the leading one, while the decay H + — > W + j induced at one-loop level has 
a BR relatively large of order 10~ 2 . However for 180 GeV< m H +, the dominant decay of the 
charged Higgs boson is the mode tb. Then, see Fig. [31(c), for the case with tan /3 = 1 one gets that 
BR(F+ -)• t + u t ) « 1 when m H + < 175 GeV, while in the range 175 GeV < m H + < 180 GeV the 
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mode ts is relevant and the BR{H + — > W + ^) even is of order 10" 4 . For 180 GeV< m H+ < 600 
GeV, the decay channel W + hP becomes relevant, whereas for the range 600 GeV < m H + the mode 
VF+A is dominant. It is convenient to mention that this sub-scenario is special for the mode tb, 
because its decay width is zero at the tree-level, since the CKM contribution is canceled exactly 
with the terms of the four-zero texture implemented for the Yukawa coupling of the 2HDM-III. 
For tan f3 = 10, we show in plot Fig. EJd) that the dominant decay of the charged Higgs boson is 
the mode t + v t , when m H + < 180 GeV, but for 180 GeV< m H + < 250 GeV, the decay channel 
tb becomes the leading one, whereas for the range 250 GeV < m H +, the mode W + h° is again 
dominant. The decay H + — > W +, y has a BR of order 10 -5 to 10~ 7 , which could be irrelevant 
channel decay for LHC. 

In order to cover further the Higgs sector in our analysis, it is appropriate to also mention how 
the previous results change with m^o, m^o and a. In general the behavior of the decay modes 
of the charged Higgs boson is similar to the cases presented above, except for the decay channel 
WhP, which can be important or unimportant. For a = 0, this mode has BR < 10 -3 when tan f3 is 
large. However, for tan/3 < 1, it becomes the dominant one. In the case a = /3, the decay channel 
WhP can be the dominant mode, with a BR of O(l). 

As a general lesson from this section, we notice a distinctive features of our 2HDM-III, namely 
that the decay modes W + j becomes very large, with a BR = 2 x 10 _1 or 10~ 2 for some of the 
scenarios considered, with tan (3 < 2. We also want to compare the general behavior of these decay 
modes, within 2HDM-III with the 2HDM-II case. In order to compare the 2HDM-III results with 
those in the 2HDM-II, we show in Fig. H the BR(H + W + j) m in the scenario A and B from 
2HDM-III, as well as BR (H+ -)■ W + -f) u from 2HDM-II vs. m H +, taking again tan/3 = 0.3, 0.5, 
1, 10. We observe that the mode W +, j is important when 170 GeV < < 180 GeV and for 

0.1 < tan j3 < 1, taking x*- d = —1. Thus, we find that the effect of the modified Higgs couplings 
typical of the 2HDM-III shows up clearly in the pattern of charged Higgs boson decay H + — > W + ^f, 
which can be very different from the 2HDM-II case, and thus enrich the possibilities to search for 
states at current (Tevatron) and future (LHC, ILC/CLIC) machines. 



IV. EVENTS RATES AT LHC FROM THE DECAY H+ -» W+>y 



The production of charged Higgs bosons at hadron colliders has been evaluated in 



for the Superconducting Super Collider, SSC) and more recently 
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25( (for the LHC) literature, 
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FIG. 4: The figure shows the BR(H+ -> W+j) uh BR(H+ -> W+j) u vs. m H +, taking tan ,3 = 
0.3, 0.5, 1, 10 for Xif = -1 (dots), xtf = or case 2HDM-II (solid line), xtf = 1 (dashes). 

mainly for the 2HDM-II and its SUSY realization (i.e., the MSSM). In these two scenarios, when 
kinematically allowed, the top quark decay channel t — > bH + is the dominant production 
mechanism. Instead, above the threshold for such a decay, the dominant H ± production reaction 
is gluon-gluon fusion into a 3-body final state, i.e., gg — > tbH^ 5 . Both processes depend on the 
coupling H~tb and are therefore sensitive to the modifications that arise in the 2HDM-III for this 
vertex. However, detection of the final state will depend on the charged Higgs boson decay mode, 
which could include a complicated final state, that could in turn be difficult to reconstruct. For 
these reasons, it is very important to look for other production channels, which may be easier 
to reconstruct. In this regard, the s-channel production of charged Higgs bosons, through the 
mechanism of c6-fusion, could help to make more viable the detection of several charged Higgs 



5 In fact, these two mechanisms are intimately related, see below. 



16 



boson decay channels [3d ]. 

A. Evaluation of event rates from direct production of charged Higgs bosons at the LHC 

The H^qq' vertex with large flavor mixing coupling, that arises in the 2HDM-III, enables 
the possibility of studying the production of charged Higgs boson via the s-channel production 
mechanism, cb — > H + + c.c. This process was discussed first in Ref. [3^] and recently by [yj] in 
the context of the 2HDM-III. 

To illustrate the type of charged Higgs signatures that have the potential to be detectable at the 
LHC in the 2HDM-III, we show in the Table [J the event rates of charged Higgs boson through the 
channel cb — > H + + c.c, alongside the corresponding production cross sections (ex's) and relevant 
BR(// + — > W + j), for a combination of masses, tan/3 and specific 2HDM-III parameters amongst 
those used in the previous sections (assuming m^o = 120 GeV, m^o = 300 GeV and the mixing 
angle at a = tt/2 throughout). In particular, we focus on those cases where the charged Higgs 
boson mass is above the threshold for t —> bH + , for two reasons. On the one hand, the scope of the 
LHC in accessing t — > bH + decays has been established in a rather model independent way. On 
the other hand, we have dealt at length with the corresponding BRs in section III. (As default, we 
also assume an integrated luminosity of 10 5 pb _1 .) In all cases we use the results of our previously 
work [141 ]. 

To illustrate these results, let us comment on one case within the scenario B, because is the 
case most conservative. From Table [U we can see that for Scenario B, with (xfj = 1jx£- = 1) 
and tan f3 = 0.3, we have that the is heavier than rat — n^b, as we take a mass m H + = 200 
GeV, thus precluding top decay contributions, so that in this case a(pp — > H + X) ~ 2.2 x 10 -2 
pb, and including the decay H + — > W + *y , in the final state we set 42 events. Even for charged 
Higgs masses as long as 300 GeV, we can still obtain 4 events. The most promising rate for decay 
H + — > W +/ y, is for mu+ = 200 and tan /? = 1, when we get 63 events. 

B. Indirect production of charged Higgs bosons at the LHC 

We have found that, in some of the 2HDM-III scenarios envisaged here, light charged Higgs 
bosons could exist that have not been excluded by current experimental bounds, chiefly from LEP2 
and Tevatron. Their discovery potential should therefore be studied in view of the upcoming LHC 
and we shall then turn our attention now to presenting the corresponding hadro-production cross 
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TABLE I: Summary of LHC event rates for some parameter combinations within Scenarios B with an 
integrated luminosity of 10 5 pb -1 , for the signal H + — > W + j, through the channel cb — > H + + c.c. 





tan (5 


m[l+ in GeV 


cr(pp — > H+ + X) in pb 


BR(H+ -5- W+j) 


Nr. Events 




0.3 


200 


2.1 x 10 2 


2 x 10~ 6 


42 


(i,i) 


0.3 


300 


3.1 x 10 


1.1 x 10" 6 


4 


(i,i) 


1 


200 


3 


2.1 x 10" 4 


63 


(i,i) 


1 


300 


0.98 


4.2 x 10" 7 







10 


200 


1.1 


9.8 x 10" 5 


10 


(i,i) 


10 


300 


5.2 x lO" 1 


1.1 x 10" 6 






sections via an indirect channel, i.e., other than as secondary products in (anti)top quark decays 
and via c6-fusion, considered previously. 

In previously work we evaluate the qq,gg — > tbH~ + c.c. cross sections. We use HERWIG 
version 6.510 in default configuration, by onsetting the subprocess IPRDC = 3839, wherein we 
have overwritten the default MSSM/2HDM couplings and masses with those pertaining to the 
2HDM-III: see Eqs. 

Now, let discuss again the scenario B. From Table HU we can see that for this scenario, with 
(xtj = 1) xfj = 1) an d tan f3 = 0.3, we have that the is heavier than nit — m bi as we take a mass 
rriff+ = 200 GeV, thus precluding top decay contributions, so that in this case a(pp —> tbH~) ~ 25.8 
pb, while the signal H + — > W +r y give a number of events of 5. When the mass of the charged 
Higgs take the value of 300 GeV, we can obtain only 1 event. One can see that for m#+ = 200 
and tan/3 = 1, we get 48 events. However, the case most promising rate for decay H + -4 VF + 7 is 
when niH+ = 200 and tan (3 = 10, we get the spectacular 122 events, which could be important for 
the LHC collider experiments. 

Altogether, by comparing the q q, g o — > tbH~ + c.c. cross sections herein with, e.g., those of the 

n 



MSSM in [6j or the 2HDM in [4J, |45j], it is clear that the 2HDM-III rates can be very large and 
thus the discovery potential in ATLAS and CMS can be substantial, particularly for a very light 
H^, which may pertain to our 2HDM-III but not the MSSM or 2HDM-II. However, it is only by 
combining the production rates of this section with the decay ones of the previous ones that actual 
event numbers at the LHC can be predicted. 



V. CONCLUSIONS 



We have discussed the implications of assuming a four-zero Yukawa for the properties of the 
charged Higgs boson, within the context of a 2HDM-III. In particular, we have presented a detailed 
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TABLE II: Summary of LHC event rates for some parameter combinations within Scenarios B with for an 
integrated luminosity of 10 5 pb -1 , for H + — > W + -f signature, through the channel qq,gg — )■ tbH + + c.c. 





tan f) 


m H+ in GeV 


cr(pp — > H+tb) in pb 


BR(H+ -> W+j) 


Nr. Events 




0.3 


200 


25.8 


2 x 10~ 6 


5 


(i ; i) 


0.3 


300 


10.1 


1.1 x 10" 6 


1 


(i,i) 


1 


200 


2.3 


2.1 x 10~ 4 


48 


(i,i) 


1 


300 


0.79 


4.2 x 10~ 7 





(i,i) 


10 


200 


12.5 


9.8 x 10~ 6 


122 


(i,i) 


10 


300 


0.48 


1.1 x lO" 6 






discussion of the charged Higgs boson couplings to heavy fermions and the resulting implications 
for the decay H + — > W + j, induced at one-loop level. The latter clearly reflect the different 
coupling structure of the 2HDM-III, e.g., with respect to the 2HDM-II, so that one has at disposal 
more possibilities to search for states at current and future colliders, enabling us to distinguish 
between different models of EWSB. We have then concentrated our analysis to the case of the 
LHC and showed that the production rates of charged Higgs bosons at the LHC is sensitive to 
the modifications of the Higgs boson couplings. We have employed the results of the s-channel 
production of through effusion and the multibody final state induced by (/(/-fusion and qq- 
annihilation. Finally, we have determined the number of events for the most promising LHC 
signatures for the decay H + — > W + ^f within 2HDM-III, for both cb —> H + + c.c. and qq — > ibH + 
+ c.c. scatterings (the latter affording larger rates than the former). Armed with these results, 
we are now in a position to carry out a detailed study of signal and background rates, in order to 
determine the precise detectability level of each signature. However, this is beyond the scope of 
present work and will be the subject of a future publication. 
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